The parametric resonance in the a^-dynamo models 



Maxim Reshetnyak 

Institute of the Physics of the Earth RAS 
B.Gruzinskaya 10, Moscow, 123995, Russia 
email: m.reshetnyak@gmail.com 

Abstract 

It is shown, that retardation in the a-quenching in the Parker's dynamo 
model leads to parametric resonance. This result is observed in the numer- 
ical simulations and can be reproduced in the simple analytic model. The 
other interesting effect in the model with retardation is a appearance of the 
long-term processes with the period much larger than the typical time of 
retardation. 
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1 Introduction 



The main idea of the dynamo theory, which is believed can explain existence of 
the magnetic fields observed in cosmos, is that kinetic energy of the conductive 
motions is transformed into the energy of the magnetic field. Magnetic field gen- 
eration is the threshold phenomenon: it starts when magnetic Reynolds number 
R m reaches its critical value R^. After that magnetic field grows exponentially up 
to the moment, when it already can feed back on the flow. Description of transi- 
tion from the linear regime when influence of the magnetic field onto the flow is 
negligible to the nonlinear regime, when magnetic energy can exceed the kinetic 
energy orders of magnitude (like it is in the planetary cores) is a subject of the 
modern researches iBrandenburg & Subramanian I (|2005|) . 

As a result, even after quenching the saturated velocity field is still large 
enough, so that R m » R^. Moreover, velocity field taken from the nonlinear prob- 
lem (when the exponential growth of the magnetic field stopped) can still generate 
exponentially growing magnetic field providing that the fee d back of the mag 



netic field on the flow i s omitted (kinematic dynamo regime ) Cattaneo & Tobias 




Tilgner (2008); Tilgner & B randenburg I ( 2008): Schrinner. Schmitt. Cameron 



Hejda & Reshetnyak I (|2010l) . It appears, that problem of stability of the 



full dynamo equations including induction equation, the Navier-Stokes equation 
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with the Lorentz force differs from the stability problem of the single induction 
equation with the given saturated velocity field taken from the full dynamo solu- 
tion: stability of the first problem does not provide stability of the second one. 
Moreover, some regimes close to the Case 1 from the geodynamo benchmark 



Christensen et al. I (1200 1|) are stable in contrast to the solutions with the periodi- 
cal bound ary cond i tions i n space and influence of the boundary conditions can be 
important iTilgner I (|2008|) : ISchrinner. Schmitt. Cameron I (|2010|) . 

One of the si mple explanations of this phenomenon (at least for some regimes) 



was offered in iReshetnyak I (|2010|) . Using Parker's dynamo model for the thin 
disk it was shown, that the ^-effect, taken from the nonlinear oscillating saturated 
problem, can still generate exponentially growing magnetic field. The origin of 
this effect is closely related to the parametric resonance. Here we develop these 
ideas and show how the phase shift in the a-quenching can effect on the behaviour 
of the generated magnetic field. 



2 Parker's dynamo 



One of the simplest dynamo models used in t he sol ar and galactic applications is 
the Parker's one-di mensional model IParker I (1 197 lb (see its development for the 
galactic dynamo in lRuzmaikin. Shukurov & Sokoloff I ([1988)): 



dA 

— = aB + A 

ot 



dB 

— = -DA' + B' 
ot 



(1) 



where A and B are azimuthal components of the vector potential and magnetic 
field, a(z) is a kinetic helicity, D is a dynamo number, which is a product of 
the amplitudes of the a- and ^-effects and primes denote derivatives with respect 
to a coordinate. For the Galaxy the only one left coordinate is cylindrical polar 
coordinate z. For the thin shells, which we will have in mind in this paper, is a 
latitude ■&. Equation (Q~|) is solved in the interval -90° < ft < 90° with the boundary 
conditions 5 = and A' = at ■& = ±90°. System (QQ) has growing solution, when 
|£)| > |2) cr |. Putting nonlinearity of the form 



<*o(fl) 

1 "h Em 



(2) 



in (Q]), where E m = (B 2 + A' 2 )/2 is a magnetic energy, gives quasi- stat ionary solu- 
tions for the positive D, see about various forms of nonlinearities in iBeck et al. 
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(|1996|) . The property of the nonlinear solution is mostly predetermined by the 
form of its first eigenfunc tion. 



The main result of iReshetnyak I (|2010f) was, that simultaneous solution of 
equations for A, B, a CD © and similar equations for the new magnetic field 
(A, B): 



dA ~ ~ 
— = aB + A 
dt 



a d 

— = -DA' + B' 
dt 



(3) 



with the same a can lead to the exponentially growing (A, B), when the field 
(A, B) is already saturated. It happens when (A, B) oscillates, and initial condi- 
tions for (A, B) and (A, B) are slightly different. T his effect is very similar t o 
what was observed in the more sophisticated models Cattaneo & Tobias I (|2009i) . 
Tilgner & Brandenburg I (|2008|) . lHejda & Reshetnyak I (|2010l) . Below we consider 
this effect in more details. 



3 Parker's dynamo with retardation 

Here we return to the system (CO)- So as a is a function of the magnetic field which 
oscillates, in the general case we are in a position to expect appearance of the 
parametric resonance for (A, B), as well. Now we consider the more general form 
of the nonlinearity with retardation r: 

a(» t r) = ^ (4) 

and find how r effects on the magnetic field generation. Results for a = sin(2#) 
and D = 300 shown in Fig. Q] are quite unexpected: in spite of the fact, that a de- 
pends on the squared magnetic field, the mean magnetic energy is not symmetric 
on r = 0.5 (r is in units of process's half period Tq « 0.45). After some decrease 
of the magnetic field amplitude for < r < r min = 0.17 magnetic field starts to 
increase being periodical, see Fig. [2] The sharp decrease at r = r br « 0.5 (see 
Fig. takes place up to the values at T min , accompanied by the long-term mod- 
ulation with the period 8 times larger than the period of the original oscillation 
of E m Tq at r = 0. The amplitude of this oscillation increases with increase of t, 
see Fig. [2]l-e. The amplitude of the new oscillation starts to change at r = 0.89, 
resembling beating, see Fig [2^. 

Increase of t up to 0. 1 leads to the change of the form of the radial B r = -A' & 
and azimuthal B components from sinusoidal to the saw-shaped form, see Fig. [3] 
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Figure 1 : Dependence of the mean over the volume and time squared magnetic 
field B 2 on retardation r in a-quenching. B 2 is normalized in such a way, that 
B 2 \ r=0 = 1. t is measured in units of T . 



accompanied with delay of B r from B. For the larger r the phase shift (p decreases. 
Increase of r to 0.33 leads to the growth of the magnetic field and appearance of 
the pike-like extremums. This region of r corresponds to the parametric reso- 
nance. 
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Figure 2: Evolution of the mean over the volume magnetic energy for the different 
retardation: r = (a), r = 0.33 (b), r = 0.57 (c), r = 0.67 (d), r = 0.89 (e). 



The further increase of r leads to the sharp decrease of the magnetic field 
amplitude, see Figs. |T1 El [3j The form of the curve is close to that one for r = 0.1 
with one exception: the new large period appears. Note, that in Yoshimura (1978a, 
1978b) such long-term modulation was riched at r > T , see review of some other 



sources of the long-term variations in the solar dynamo in iTobias I (120021) . The 
different behaviour of E m for < r < 0.5 and 0.5 < r < 1 may be explained 
by the accumulation of the time shift between the magnetic field and a. This 
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Figure 3: Evolution of the mean over the half of the volume B r (left column) and 
B (right column) for r = - a, b; 0.1 - c,d; 0.33 - e,f; 0.57 - g,h; 0.67 - i,j; 0.89 
-k,l. 



suggestion is supported by the fact that for r > 1 behaviour of E m is very close to 
that one for r > 0.7. 

Up to the moment we did not consider how the magnetic field depends on 
§. The dipole solution is a wave propagating from the poles to the equator with 
maximum at the middle latitudes for t = 0, see Fig. |4j Increase of t leads to 
the stripe-like solution, i.e. more contrast changes of the sign of the fields and 
increase of the magnetic field at the poles Fig. @}). The long term variations are 
well resolved in Fig. HI:. 
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Figure 4: The butterfly diagram B(t, 6): r = - a; 0.33 - b and 0.67 - c. 

4 Parametric resonance 



To consider how parametric resonance appears we follow iReshetnyak I (|2010f) and 
present solution for (A, B) in the form of the waves: B = bs'm(t),A = sin(t+(f+9), 

B = sin(? + 6), and a = we get how generation depends on r. Then 

1 + B 2 (t - t) 

putting it in (QQ) we get two equations. Equation for B does not include r, so we 

In 

consider only production of A 2 . Then 6A((p, 6) = ckq J ^ 2 dt. If |II| » 1, 

o 1 + 

6A(w, 9) 

where n = , then (A, B) is unstable. 

6A((p, 0) 

The integral for 8 A gives: 

6A((p, T) = hi+ h 2 ig(<p), h Y *\- 0.3 cos(r) 2 , 
Then 



h 2 ss 0.8 sin(2r). 



n = i + 



0.8 sin(2r) 



rtg(^) 



(5) 
(6) 



1 - 0.3 cos(r) 2 
n 

and for cp — > ±— (what corresponds to the phase shift between the components 



fB r B ~ 0) and t ^ |II| grows, and parametric resonance appears, see Fig. [5J 
Note, that this analysis explains small decrease of B 2 at small r (see Fig. Q]) which 
corresponds to the positive values of tp in Fig. [5] 
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Figure 5: Dependence of II(r) for different (p: - solid line, -85° - circles, -70° 
- triangles, 70° - squares, 85° - stars. The straight line corresponds to II = 1. 



5 Conclusions 

Here we considered the only one form of the nonlinearity for the fixed value of the 
dynamo number D as a function of the time lag r in the a-quenching. However 
even this simple model demonstrates variety of effects: change of the form of the 
poloidal and toroidal fields, regions of the weaker and stronger fields, appearance 
of the long-term variations. The quite natural choice of r can lead to the sharp in- 
crease of the magnetic field amplitude concerned with the parametric resonance. 
This explanation does not contradict to the simple linear analysis presented above. 
It is very tempting to correspond appeared in simulations the long-term periodic- 
ities with that ones of the solar activity larger than the main 22 years period. The 
difficulty is to justify the choice of the particular r which should be obtained from 
the solution of the more sophisticated nonlinear model. 




I thank D.Sokoloff for discussions. 
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